Ionic phase transitions in non-ideal systems 
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We construct an explicitly solvable Landau mean-field theory for volume phase transitions of 
confined or fixed ions driven by relative concentrations of divalent and monovalent counterions. 
Such phase transitions have been widely studied in ionic gels, where the mechanism relies on self- 
attraction or elasticity of a network. We find here that non-ideal behavior of ions in aqueous solution 
can in theory drive phase transitions without a self-attracting or elastic network. We represent non- 
ideality by a Debye-Huckel-like power-law activity, or correlation free energy, and retain a mechanical 
self-repulsion to avoid runaway collapse due to the non-ideal term. Within this model we find a 
continuous line of gas-liquid- type critical points, connecting a purely monovalent, divalent-sensitive 
critical point at one extreme with a divalent, monovalent-sensitive critical point at the other. An 
alternative representation of the Landau functional handles the second case. We include a formula 
for electrical potential, which may be a convenient proxy for critically varying volume. Our relatively 
simple mean-field formulation may facilitate explorations of tunable critical sensitivity in areas such 
as ion detection technology and biological osmotic control. 

PACS numbers: 82.60.-s, 05.70.Jk, 87.15. Zg, 64.70. Nd, 82.60.Fa, 64.70.fh 



In charged polymer networks, phase transitions were 
observed and explained theoretically some time ago by 
Tanaka and co-workers [IH1]. In experiments, gels can 
exhibit ionically-driven collapse or expansion as the os- 
motic pressure responds to slight changes in external ion 
or composition. Volume phase transitions have long been 
discussed with an eye to applications [5J |B] and as candi- 
date mechanisms for essential biophysical processes [TJ. 
In previous models network self-attraction and restoring 
elasticity has played a crucial role, invoking specific elas- 
tic properties of polymer networks [8] ■ 

Here we seek to recast the theory of ionic phase tran- 
sitions into a simple and mathematically solvable mean- 
field formulation, including (in simplified form) the non- 
ideality that is present in all ionic solutions. We find 
that in principle a self-attracting network is not neces- 
sary for a discontinuous phase transition in the presence 
of a simplified Debye-Hiickel-type non-ideality, because 
non-ideality itself acts as an effective self-attraction that 
mediates an ionically driven phase transition of the gas- 
liquid type. Our model requires us, however, to attribute 
a self-repulsion to our fixed-ion system to avoid runaway 
collapse due to the non-ideal term. 

Our simplified formulation of ionic transitions may, in 
practical terms, provide a tunable critical sensitivity to 
ion valence or concentration. This suggests, for example, 
a mechanism for biophysical cellular functions such as 
homeostasis. The theory may also lend itself to engineer- 
ing applications involving ion detection. For these pur- 
poses electric potential <&, whose behavior is predicted, 
might often be a more convenient dependent parameter 
than volume in applications to ion detection or in bio- 
physical roles for critical ionic sensitivity. 

Consider a population of N ions each of charge go 
that may be bound to a mechanical structure, such as 
a polymer network, or otherwise confined within a vol- 



ume permeable to counterions. These are the conditions 
for Donnan equilibrium [5j in which osmotic pressure of 
excess counterions goes hand-in-hand with an internal 
voltage relative to outside, with the same sign as go- 
Even for $ ^ 0, neutrality holds to a good approximation 
whenever the voltage drop $ occurs only at the bound- 
ary, or more generally if N 3> C$/g , where C is the 
system capacitance. 

The immobile ions, confined within a variable volume 
V, have concentration c = Nq/V — 1/v. Monovalent 
and divalent counterions with charges q a = —go and 
g5 = — 2g are introduced at external concentrations a 
and b, and have within V the concentrations a' = N a /V 
and b' — Nb/V. Their free diffusion in and out is con- 
trolled by the the electrical potential We introduce 
nonidcality of mobile ions by way of a Debye-Hiickel in- 
teraction or correlation free energy |10| 

-Fdh = -v^k B TV(j:zlc a ) 3/2 (1) 

a 

where q a — z a e. The essential features of -F DH are that it 
is negative and contains mobile-ion concentrations raised 
to the 3/2 power. The constant v B is proportional to the 
cube of Bjerrum length. In standard aqueous conditions, 

— = (12vr) 2 (^Y « 3.4 molar. (2) 
v B V e z ) 

To obtain closed-form solutions, we omit confined ions 
from the sum, and replace the factor (a' + 46') 3 / 2 by (a' + 
26') 3 / 2 , which will equal c 3 / 2 once neutrality is imposed. 
Thus we employ the Debye-Huckel-like interaction term 

F iut = -v^k B TV(a' + 2&') 3/2 . (3) 
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With a' — N a /V and b' = N b /V and taking derivatives, 
-dF int /dV = P int = -i< 2 k B T (a' + 2b') 3 / 2 (4) 
dF int /dN a = k B T\n la = -§ k B Tv B /2 (a' + 2b') 1 ' 2 (5) 
dF int /dN b = k B T\n lb = -3v B /2 k B T (a' + 2b') 1 / 2 (6) 



Here 7 a and 7b are activity coefficients [TT] for the mono- 
valent and divalent ions: /i a = [J,® + k B T\nj a a. The 
proportionality of In 7 to the square root of ionic con- 
centration represents Debye-Hiickel behavior, a theory 
appropriate to low ionic strength |12) . The free energy, 
with $ externally controlled, is 



F(N a ,N b ,V,$) 
+ F Q (V) 



F iDt (N a ,N b , V) + (N -N a -2N b )q $ 

Nn . N b 



k B T 



N a \n 



c eV 



N b ln- 



c eV 



(7) 



(with e the natural log base). Here Fo(V) describes me- 
chanical constraints on the immobile ions, such as a poly- 
mer network carrying the fixed ions or a membrane con- 
taining them. The logarithmic terms are the ideal free 
energy of the mobile ions, and the concentration scale 
Co will cancel out. The mechanical contribution to pres- 
sure is — dFo/dV = Pq (equivalently a function of v or 
c) and the overall pressure and the chemical potentials 
[i a = dF/dN a and fi b — dF/dN b are 

P = P Q (v) + (a' + b')k B T + P int (8) 
Ma = -go* + fc B T , ln(7 Q a / /co) (9) 
M 6 = -2g $ + fe B Tln(7 6 &7co) (10) 

If we set [i a = fc B Tln(o/co) and /j, b = k B T ln(6 / cq) , a 
and b become effective concentrations, referred to an 
ideal external solution. To fix $, we impose neutral- 
ity, c — a' + 2b' , and introduce a dimensionless potential 
(f> = qo$/k B T which will always be positive (since <!> will 
always be of the same sign as qo). We have altogether 

P = P (N /c) + \{a! + c)k B T - \vT <?< 2 k B T (11) 

= ln(o7o) - c 1 ' 2 (12) 

20 = ln(67&) - H /2 c 1/2 (13) 

which, together with c = a' + 2b' = Nq/v, we will solve to 
find P[v). Combining the equations involving 4>, putting 
2b' + a' — c = for neutrality, and solving yields 

a' = {a 2 /4b)[(l + (8b/a 2 )c) 1/2 - l] (14) 
We define the dimensionless divalent parameter 



(3 



(15) 



To motivate use of /3, imagine n divalent ions B coopera- 
tively exchanging with 2n monovalent ions A, according 
to CA2 n + nB ^ CB n + 2nA. Equilibrium is then 



[CB n ] 
[CA 2n ] 



= K 



A] 2n 



(16) 
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FIG. 1: Coexistence diagram in dimensionless volume and 
pressure for a system tuned by repulsion a to the monovalent 
(A = 1) critical point. The divalent parameter ft rises from 
zero upon addition of divalent ions, so that j3 < is not 
possible. In similar coexistence diagrams for < A < 1, the 
upper region will be accessible. 



For large n, the discontinuous transition CB n —¥ CAm 
occurs at a fixed value of b/a 2 . 

The dimensionless pressure p(v) = Pv B /k B T is 



Po{v) + 



1 v B 
27 



1/2 



- 1 



1 /W B \ 3 / 2 

2 VV 



(17) 



where po = Pov B /k B T. If the system is to be stable 
against collapse to v = due to the interaction term, 
we require that po(v) include a repulsion diverging faster 
than l/i! 3//2 . We will use the minimal choice 



Po(v) = +a(v B /vf 



(18) 



(with a dimensionless). The dimensionless electrical po- 
tential is 



= In 



1+13- 



1/2 
- 1 



3 /v R \ 1 /i 



ln{v B a). (19) 



From the Gibbs-Duhcm relation (at constant T, or here 
/3) we obtain the chemical potential, linking coexisting v 
and v' at specified /?, as fx(v) — pv — J p(v)dv, or 



H(v) = 2ajv-\v 1/2 -In v+y/v(v+/3) 



(20) 



Self-intersections of the curve (p(v),/i(v)) yield phase 
boundaries as in Figs 1 and 2. 
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FIG. 2: Critical line (dotted) connecting the monovalent 
(A = 1) and divalent (A = 0) critical points. The divalent 
parameter is (3 and p is dimensionless pressure. Selected crit- 
ical points and their phase boundaries (solid) are shown. The 
phase boundary becomes inaccessible (moving to (3 > oo) for 
the critical point at A = 0. 



We solve our model in terms of the parameter 

x = \/v B /v. (21) 

At a critical point the conditions p'(v) = p"(v) = and 
p'(x) = p"{x) = are equivalent. Regarding the critical 
value of p as a function of x and /3, and introducing the 



parameter A = /3 c x^, , we have the ^-derivatives 

p x = x c [4ax 2 c - §a; c + l +A] = (22) 
p xx = 12a x\ - 3x c + 1 + A 3 = (23) 

Eliminating a in favor of A gives the line of critical points 

„2 r i 



[^(1-6A + A 3 ) + A/(1 + A)] (24) 

(25) 
(26) 
(27) 



x c = |[2 + 3A-A 3 ] 
a c = (1 + 2A- A 3 )/4x 2 
(3 C = 1/(.^A 2 ) 



The interval < A < 1 corresponds to oo > /3 > 0. Al- 
though we view a(A) as the mechanism by which ionic 
critical points can be tuned, we can also view a c — a(A) 
as a critical value of the repulsion strength. With ionic 
conditions fixed, we could drive transitions by modulat- 
ing a, analogous to transitions in ionic gels induced by 
variation of solvent composition [5] . 

About each critical point, with Sv = v — v c and 8(3 = 
(3 — (3 C , we construct a Landau expansion for the pressure, 



p((3,e) = Pc -AS/3 + BS/3e 



Ce 3 



(28) 



where e — 8v + n8f3 is an order parameter linear in Sv 
and 8(3. In a potential application, we imagine tuning a 
system to a critical point, (3 — f3 c and p — p c . Within 
mean-field theory, when the divalent ratio is changed by 
8(3 = (3 — (3 C one predicts a singular expansion or contrac- 
tion 



-{A/C) 1/3 (<5/3) 1/3 



(29) 



where the cube root is taken with the same sign as 8(3. As 
discussed below, this singular behavior in volume might 
in practice be better monitored via the electric potential 
than the volume change Sv or e. 

To evaluate the Landau coeffients, we expand the pres- 
sure around the critical point, to third order in Sv = v—v c 
and to first order in 8(3, giving 

p{v,(3) ss p c + pp8(3 + p v p8(3e + \p vvv e 3 (30) 

where (3 and v subscripts denote partial derivatives and p, 
p v , etc. are critical values. We have chosen k = p vv /3 /p vvv 
to eliminate an e 2 term. Evaluating the various deriva- 
tives, the Landau parameters along the critical line are 



A = 
C = 



±^A 3 /(1 + A) 2 , B = \xl\ 3 
^a;8[2 + 12A-7A 3 + 3A 5 ] 
(A 3 /x c )(A 2 + l)/(8a; c a c -/? c A 5 ) 



(31) 
(32) 
(33) 



where x c and a c are known functions of A from Eqs ( 25 ) 



and (26 1. C has no zeros within < A < 1, while k 
is zero only at A = 0. At the monovalent critical point 
{(3 = 0, i.e. A = 1) Eqs ([24} through ([27} yield p c = 
32/27 « 1.185, x c = 8/3 w 2.67, a c = 9/128 0.0703, 
the Landau coefficients are 



.4 
B 



2 9 /3 4 ■ 
2 16 /3 6 



6.321 (34) 

» 89.90 (35) 

C = 5(2 21 /3 8 ) « 1598.2 (36) 

and the order parameter is e = Sv + (1/2)8(3. Here 
the system enters coexistence for any value of /3 > 0, 
with e R3 —(0.158) (3 1 / 3 . A system tuned to the monova- 
lent critical point moves with infinite response towards 
a smaller volume with the introduction of any divalents. 
In the liquid-gas analogy, this path follows the critical 
isobar. 

At the purely divalent critical point ((3 = oo, A = 0) 
both A and B vanish, but we can here characterize the 
line of critical points by the alternative parameterization 



p((3,e) = p c + A>8(p-^)-B>8(l3-^) 
A' = z c (l-A) 3 / 2 /(l + A) 2 
B' = \x 3 c (l-\ 2 ) 3 / 2 



Ct 



(37) 
(38) 
(39) 



The parameter (3 1/2 is proportional to the monovalent 
concentration at fixed divalent concentration. The vol- 



ume singularity ( 29 ) can be rewritten as 
-A' 



C 



1/3 



(40) 



4 



where again the cube root has the same sign as its argu- 
ment. At the divalent critical point (/J = oo, A = 0) Eqs 
([24]) through (|7| yield Pc = 4/27 « 0.1481, x c = 4/3, 
a c = 9/64 w 0.1406, 
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A' = 4/3, £?' = 32/27 
C = 2 13 /3 9 « 0.4162 



1.185 



(41) 
(42) 



For this divalent critical point at /3,T 1/2 = 0, the order pa- 
rameter is simply the volume, e = 5v. A purely divalent 
system that is tuned to be critical moves with infinite 
response towards a larger volume with the introduction 
of any monovalents, as Sv ss (1.214) /3 -1 / 6 . 

Returning to the dimensionless potential <j) = qo$/k B T 



in Eq ( 19 ), along the critical line with a constant we find 



dv 



3z£ 
4t> B 



1 + 2A - X 6 

2 + 3A- A 3 



(43) 



Since S(f> ~ <t>v$ v i the potential will display the same 
power-law singularity as the volume. 

Our simplified framework is intended to motivate in- 
vestigation of ionic phase transitions in ever simpler sys- 
tems. As a fascinating analogue of liquid-gas transi- 
tions, and in the context of ionic gels, such transitions 
have been a topic of discussion for many years. Possi- 
ble roles in biological systems [7] and potential applica- 
tions such as metal ion detection and electrical sensitiv- 
ity [3 [6] should remain applicable to phase transitions as 
discussed here. 

The authors would like to thank colleagues at the De- 
partment of Physics and Astronomy at Washington State 
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